Early reionization by cosmic strings revisited 
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Measurements of the CMB temperature anisotropy and the temperature-polarization cross cor- 
relation by WMAP suggest a reionization redshift of z ~ 17 ± 5. On the other hand, observations 
of high redshift galaxies indicate a presence of a significant fraction of neutral hydrogen at redshift 
z ~ 6 — 7. We show that cosmic strings with tensions well within, but not far from, current obser- 
vation bounds could cause early star formation at a level sufficient to explain the high reionization 
redshift. 
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Most of the baryonic matter in the present day uni- 
verse is in the form of ionized plasma. Yet we know that 
the universe was neutral immediately after the recombi- 
nation. Therefore, sometime between now and redshift 
z ~ 10 3 the universe was reionized. The detailed history 
of reionization is a mystery and has become a particularly 
hot topic of discussion following the release of WMAP's 
first year data. The WMAP team has reported an opti- 
cal depth to the last scattering r = 0.17 ±0.06, which for 
the instantaneous and homogeneous reionization implies 
a redshift of z ~ 17 ± 5 0, 0- The precise value of the 
reionization redshift inferred from CMB data depends on 
the assumed ionization history and a more conservative 
interpretation of WMAP's result is that the measured 
optical depth to last scattering is consistent with a sig- 
nificant ionized fraction at redshifts z ~ 11 — 30 Q{. 

The commonly accepted general picture of reionization 
is that it was caused by first stars and quasars emitting 
high energy photons. In order for the universe to be 
significantly ionized by redshift z r ~ 20 a sufficiently 
large fraction of matter would need to be collapsed to 
form stars prior to z r . It remains to be seen if collapse 
of Gaussian initial inhomogeneities can produce the re- 
quired number of reionizing stars formed by such early 
time 0. 

To add to the puzzle, the WMAP value of the re- 
ionization redshift is significantly higher than the value 
of z r ~ 6 — 7 which one would infer from measurements 
of Lyman-a absorption in the spectra of high redshift 
galaxies (e.g. see Q). This could mean that the ionized 
fraction evolved non-monotonically with redshift and in- 
volved processes more complex than previously thought 

sauna. 

In this paper we investigate the possibility that early 
reionization was triggered by stars formed in the wakes of 
moving cosmic strings. The idea that strings could cause 
reionization has been discussed in a number of papers 
[m Il2l Il3| ] . Refs. 0,0] were written over a decade ago 
and used cosmological parameters that were significantly 
different from the ones we know now. Also, at the time 
those papers were written, cosmic strings were considered 
a prime candidate to explain the formation of large scale 
structure and were generally attributed a higher value of 



mass per unit leng th then currently allowed by observa- 
tions. In Ref. [13] stings are represented by a spectrum 
of linear perturbations. It is not clear to us why this 
should give an adequate description, since we know that 
the inhomogeneities induced by strings on scales relevant 
to reionization are highly non- linear in character . Here, 
we re-investigate the reionization by strings in the light 
of the current data, and with nonlinear string dynamics 
fully taken into account. 

Another motivation for this study is the recent real- 
ization that cosmic strings (or D-branes) are naturally 
produced in many brane inflation scenarios, motivated 
by String/M Theory [H [HE In these the- 
ories strings can move and interact in extra dimensions, 
in addition to the observed three spatial dimensions. In 
particular, while appearing to intersect in our three di- 
mensions, they may actually miss each other in the extra 
dimension(s). Hence, the effective intercommutation rate 
of these strings will generally be lower than one. As a 
consequence, one would expect more strings per horizon 
in these theories, with more small scale structure accu- 
mulating on the strings 18]. This can have interesting 
observational consequences which must be investigated. 

Let us now consider the problem in more quantitative 
detail. The approximate fraction of total baryonic mat- 
ter that needs to be in stars in order to ionize the rest 
of the gas can be roughly estimated by accounting for 
the following main effects. The ionization energy of hy- 
drogen is 13.6 eV, while nuclear fusion in stars produces 
~ 7 MeV per hydrogen atom. Not all photons emitted 
by stars have energies above 13.6 eV and each hydro- 
gen atom can recombine more than once. The number 
of ionizing photons produced per baryon depends on the 
mass and the composition of the star which produced 
them. Massive, metal-free stars (the so-called Popula- 
tion III start) can produce roughly 10 5 ionizing photons 
per baryon. Locally observed stars, which generically 
have a much lower mass and some metal content, can 



1 Thi s is confirmed by numerical studies of cosmic string wakes 
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produce ~ 4000 ionizing photons per baryon [6j. Only 
a small fraction of the gas in the halo will have time to 
form Population III stars, since those will soon explode as 
supernovae and prevent the formation of other supcrmas- 
sive stars. Hence, star formation regions quickly become 
dominated by low-mass stars. The fraction of matter 
that would need to be in stars in order to reionize the 
entire universe can roughly be estimated 0, |20L l2l| to be 

/stars ~ 1(T 3 - 10- 4 (77/10) , (1) 

where r\ (typically > 1 for z > 7) is the number of ionizing 
photons needed per hydrogen atom. In the lower 
number corresponds to metal free stars. However, the 
metallicity is likely to be increased rather quickly as the 
first stars explode as supernovae, hence we are going to 
use a more conservative estimate of /stars j corresponding 
to the higher number in Q. 

What if WMAP's estimates are true and the universe 
was indeed reionized around rcdshift z ~ 20? Assuming 
that cosmic structures grew by gravitational instability 
from initial Gaussian density fluctuations, one could ap- 
ply the so-called extended version of the Press-Schechter 
model [2^ and estimate the fraction of matter that col- 
lapsed and formed stars prior to z ~ 20. Some analytical 
and numerical work estimating that fraction has recently 
been done in Q . The fraction of matter in stars is related 
to the quantity Pb(M m i n ), which is the fraction of mat- 
ter in halos of mass Af m i n or greater, where M m i n is the 
minimum mass that a halo must have in order to form a 
galaxy. For M m ; n w 7 x 10 5 M Q Barkana and Loeb find 
Fa(Mxoin) ~ 4 x 10~ 5 jj], which is roughly consistent 
with the needed fraction given in eq. However, the 
precise estimates of / using this method and the use of 
Press-Schecher itself at such high redshift is still a matter 
of some debate. 

We are going to suggest that cosmic strings could pro- 
vide the fraction which is just as large and hence play a 
significant role in the early reionization. We will discuss 
current observational bounds on cosmic strings later in 
the paper, but only mention now that the strings we con- 
sider have tensions that are consistent with the most re- 
cent observations of CMB and large scale structure (LSS) 
and, therefore, do not contribute appreciably to structure 
formation on large scales. 

During the radiation and matter dominated eras the 
string network evolves according to a scaling solution 
[23L 124. |25| which on sufficiently large scales can be de- 
scribed by two length scales. The first scale, is the 
coherence length of strings, i. e. the distance beyond 
which directions along the string are uncorrelated. The 
second scale, L(t), is the average interstring separation. 
Scaling implies that both length scales grow in propor- 
tion to the horizon. Cosmic string simulations indicate 
that ~ t, while 

L(t) = yt , (2) 

with 7 « 0.8 in the matter era [U E3- The so-called 
one-scale model [2^, [23 , in which the two length scales 



are taken to be the same, has been quite successful in de- 
scribing the general properties of cosmic string networks 
inferred from numerical simulations. These simulations 
have assumed that cosmic strings would reconnect on ev- 
ery intersection. It is of interest to us, however, to also 
consider the case when the reconnection probability is 
less than one. Then, because of the straightening of wig- 
gles on subhorizon scales due to the expansion of the 
universe, one would still expect £(t) ~ t, but the string 
density would increase, therefore reducing the interstring 
distance. Hence, smaller intercommutation probabilities 
imply smaller 7. 

Numerical simulations show that long strings possess 
significant amounts of small-scale structure in the form 
of kinks and wiggles on scales much smaller than hori- 
zon. To an observer who cannot resolve this structure, 
the string will appear to be smooth, but with a larger 
effective mass per unit length ft and a smaller effective 
tension T . An unperturbed string (with \i = T) exerts 
no gravitational force on nearby particles. In contrast, a 
wiggly string with ft > T attracts particles like a massive 
rod. The effective equation of state of a wiggly strings 
is ftT = ft 2 [2j| and the velocity boost given by a 
moving wiggly string to nearby matter is |3(1 l3l| 

u % = iirGfiVsjs + y= '- , (3) 

where v s is the string velocity and 7 S = (1 — v 2 )^ 1 / 2 . 

Let us now consider a wake formed behind a string 
segment of length £(t) ~ t that travelled with a speed v s 
at some early time for a period of time ~ t{. For 
now, let us assume that t{ > t eq , where t eq is the time of 
radiation-matter equality. There would be wakes formed 
during the radiation era as well, and we will discuss them 
separately later in the paper. The length l w [z), the width 
w w (z) and the thickness d w (z) of the wake will evolve 
with redshift as [H 0111 

L w ~ ti , w w ^ v s ti , cl w ~ Uiti 1 I , (4j 
z z \ z J 

as long as 

z > ZiUi/v s . (5) 

At smaller redshifts the wake thickness d w becomes com- 
parable to the width w w and the wake takes on a shape 
of a cylinder whose diameter grows as z~ 3 / 2 . We find 
that the condition (JSJl is satisfied for all wakes that have 
a chance to play a role in the early reionization. Nu- 
merical simulations [24j show that average string veloc- 
ities on scales comparable to the horizon are of order 
v s ~ 0.15. The effective mass per unit length of these 
strings is ft « 1.6 fi, which implies that the second term 
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in eq. Q will dominate 2 : 



2TtG(ia 



(6) 



where we have defined a = (p — T)/fi. 

As the universe expands, the fraction of matter in the 
wakes grows in proportion to the scale factor. For wakes 
formed at rcdshift Zi, this fraction is 



fw Zi) 



2irGpa r y 



(7) 



where we have used Eq. (01 for the average inter-string 
distance. To make quantitative estimates easier, let us 
re- write this using some characteristic values for the pa- 
rameters: 



where we have used z eq w 3400. Whether any part of 
this fraction collapses into luminous objects depends on 
the values of Zj, z and Gfia. 

What are the conditions leading to the formation of 
luminous objects within a wake? As it grows, the wake 
will fragment (due to the wiggliness of the string that 
produced the wake and due to intersections with smaller 
wakes produced at earlier times) into chunks of size com- 
parable to the wake thickness d w 0, H|| • These frag- 
ments form the CDM halos into which the baryons will 
fall after decoupling from photons. Generally, one would 
expect the baryonic gas inside the wake to collapse when 
its Jeans length becomes smaller than the thickness of 
the wake. The Jeans length is defined as 



Lj(z) = c s (z) 



Gp bw (z) 



1/2 



(9) 



where c s {z) and pb w {z) are, respectively, the speed of 
sound and the energy density of baryons inside the wake. 
The infall velocity grows with time as [3{| 



( \~ 2 ( Zi Y /2 



(10) 



After a certain time the infall becomes supersonic and 
shocks develop on either side of the baryonic wake 3 . The 



2 Strings with lower inter-commutation probabilities are expected 
to accumulated even more small-scale structure due to the sup- 
pression of loop production llSl . 

3 One could allow for the possibility for the Jeans condition to 
be satisfied before the infall into the wake becomes supersonic. 
However, the condition for supersonic infall, c s < u(z), and the 
Jeans condition, roughly c s t < d w (z), are essentially the same 
inequality. Hence, it is sufficient to consider only the wakes that 
form shocks. 



compression of the gas within the shock will heat it to a 
temperature [37j 



Tshock « —m H u(zY 
16 



(11) 



where m# is the hydrogen mass. This implies that 
the baryonic speed of sound inside the wake, c s ~ 
■\/2T s hock/wff; is comparable to u{z): 



c s (z) w 0.6u(z) 



(12) 



The density of baryons inside the shock is enhanced by 
a factor of 4 compared to the background baryon density. 
In addition, if the gas is able to cool to some equilibrium 
temperature T coo i, its density would be enhanced by an 
additional factor T s hock/7c oi. Taking these effects into 
account, the baryon density in the wake can be written 
as 



Pbw(z) 



4Xb Tlliock f z 

6wGtf r c0 oi \Zi 



(13) 



where Xb = f^/fi^/ is the fraction of matter in baryons. 
Substituting eqns lfTU|) . lfT2l and (|T3|) into © gives 



Lj(z) ~ X, 



-1/2 / Tcooi 



1/2 



(f) ; 



v ^shock 

The collapse condition d w {z) > Lj(z) then leads to 



X, 



-1/2 



T, 



cool 



hock 



1/2 



< 1 



(14) 



(15) 



or. 



hock 



t. 



>6 



cool 



(16) 



where we have used X\, « 0.16. From eqns. 111|) . 110|) 

and it follows that 



hock 



37r 2 m/f / Gpa 
25 [~^T 



z 



The lowest temperature to which the gas can cool de- 
pends on its metal abundance. The metal free primordial 
gas can cool via atomic transitions of hydrogen and he- 
lium down to 10 4 K, while molecular hydrogen can cool 
the gas down to 200K 38]. Once the gas is enriched 
with metals, it could in principle cool down to the CMB 
temperature. Since the metallicity is likely to increase 
quickly, as discussed earlier, we are going to adopt 

T coo i(z) = T cmb (z) = 2.726(1 + z)K . (18) 

The collapse condition in eq. (|16f) can now be written as 

' >0.03(^) (^) , (19) 



\GpaJ V0.15 
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where the superscript (m) on Zj denotes matter era. This, 
in turn, implies a constraint on the value of Gjx (using 



Am) 



< Z e 



3400): 



P V20/ V0.15/ 



(20) 



Formation of gaseous objects can only occur after the 
recombination, because of the Compton drag. We are 
mainly interested in the formation at much later red- 
shifts, with ineq. i|20|) providing the necessary condition. 

So far, we have only discussed wakes formed after t eq . 
Let us now consider wakes formed at some time U during 
the radiation era. For these wakes the gravitational insta- 
bility only sets in at t ~ t eq . At earlier times, the surface 
density of the wakes decreases as i -1 / 2 , while the fraction 
of dark matter accreted onto all wakes formed within a 
Hubble time of ij remains roughly constant j3^| . The in- 
fall velocity into these wakes will decrease with time as 
(1 + z) until the matter starts to dominate. Therefore, 
for a wake formed at some ti < t eq , the infall velocity at 
t > t eq will be given by 



■•<" - 1* (^r fe) ■ < 2i) 



Consequently, the temperature inside the shocks that 
would form in such wakes at t > idee would be given 
by 

~ 3n 2 m H ( G[ia \ 2 z 3 
T shock {z,Zi) « — — — — \ (22) 



25 \ v s J z zf 

The necessary condition given by eq. 1)16(1 will now be- 
come 



55000 



20\ ( G(j,a \ ( 0.15 



(23) 



„(r) 



> z eq to re-write this as a lower 



Again, we can use z, 
bound on the value of G/i, which leads to the same con- 
straint as in eq. (|20|l . 

Ultimately, we want to find the fraction of matter in 
all wakes, formed before and after t eq , that satisfy the 
collapse condition in eq. ((16(1 . This total fraction can be 
expressed as f w (z,Zi) from eq. © evaluated at z, = z eq 
and multiplied by a factor a w . This factor can be written 
as a sum of the contributions from wakes formed before 
and after t 



■-eg- 



(24) 



Because the fraction of matter in wakes formed within 
each Hubble time (defined as the time required for the 
horizon to double in size) remains roughly constant dur- 

(r) 

ing the radiation era, a w is simply equal to the number 

of Hubble times between z eq and z, 
by inequality (|2*3l . Namely, 



,0) 



(r) 



which is bounded 




< 8 + 21og 2 



0.15 



(25) 



The fraction of matter in wakes formed after t eq decreases 
with the redshift as z~ l . Therefore, the dependence of 
aw on the number of Hubble times between z eq and 



z l - m ' ) is sufficiently weak and we can simply set a^j 



1. 



Hence, the total fraction of matter in the wakes that 
satisfy the collapse condition can be written as 



fw(z) 



10 _3 7 _2 a^ 



2Q\ (G^a\ 



where 



9 + 21og 2 



(G[ia\ /0.15 
i 10^ l~ 



(26) 



(27) 



Note that a w decreases with a decrease in Gfi or an in- 
crease in z, e. g. a w ~ 10 for G/ia ~ 10~ 6 and z ~ 20 
yet a w ~ 1 for G[ia ~ 10 -7 and z ~ 30. 

The condition in eq. ((16(1 does not say anything about 
the ability of the gas to cool. Only some fraction / s 
of the total number of baryons satisfying this condition 
will actually be able to cool and form stars. At present, 
there is no good theory of this fraction. However, one 
could assume, based on the current ratio of the average 
mass density in stars to the total baryon density, that this 
fraction is of order 10% |3|j . Finally, we can estimate the 
fraction of matter in stars, formed in the wakes behind 
cosmic strings, as 



io^V 2 ^ 



20 



A 
o.i 



/ G\xa\ 



(28) 



This, together with eq. (JTJ and the WMAP's upper 
bound on the reionization redshift (z r < 30), can be used 
to put an approximate upper bound on the value of Gfi: 



Gfi < 10- 



1 — 1 — 1 

7 a a w 
0.1 



0.1 



(29) 



which is roughly an order of magnitude higher than the 
upper bound reported in [lj| . Inequalities ((29(1 and 1(201) 
define the approximate range of values of G\x for which 
strings may play a role in the early reionization. In- 
terestingly, this range is quite narrow - only an order 
of magnitude wide. However, this range also depends on 
parameters a and 7, which could be significantly different 
from their canonical values if, e. g., st ring s were D-branes 
formed at the end of brane inflation [la, uR Oil • 
Current CMB and LSS data imply that Gfi < 10" 6 
[rH l40l l4l| 4 . This bound, however, also depends on the 
number of strings within a horizon or, more generally, 



4 These constraints only apply to local cosmic strings cons ider ed 
in this paper. For recent constraints on global strings see |42| . 
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on the scaling parameter 7. When this dependence is 
taken into account, the constraint based on CMB and 
LSS power spectra becomes 

G[i < 1(T 6 7 . (30) 

Substituting this into eq. I|28[l gives 

/w < 10 -3(l^)( ? )(A), (31) 

which is of the same order of magnitude as the fraction 
needed to reionize the universe. 

It may be worth mentioning that recently there has 
been a detection of galaxy lensing that looks very much 
like a lensing by a cosmic string It remains to be 

confirmed, but assuming it is a cosmic string, the implied 
value of the string tension would be G[i ~4x 10 -7 . For 
this value of G/x, and assuming that a ~ 1, from eq. i|27|) 
it follows that a w « 7 at z ~ 20. Recent analysis of 
the WMAP and SDSS data H2 shows that data may 
actually prefer a non-zero contribution from strings, cor- 
responding to Gfi ~ 4 x 10~ 7 7. The combination of these 
two measurements would imply 7 ~ 1, which for z ~ 20, 
would give 

f w « 3 x 10- 4 ( A^ , (32) 

which is close to what would be required for the reion- 
ization by redshift z ~ 20. 



Our analysis was based on the expectation that the 
metallicity of the gas in the wakes will quickly reach the 
level sufficient to prevent the formation of heavy metal- 
free stars (the so-called Population III stars). If, however, 
a significant number of Population III stars managed to 
form prior to metal enrichment, then our results would 
be modified. In particular, a smaller fraction of matter in 
the wakes would satisfy the collapse condition l|ltj|> due 
to the higher value of T coo \ for a metal free gas. On the 
other hand, a smaller value of f rm star would be required 
to reionize the universe. In combination, this is likely to 
reduce the upper bound on G/j, (eq. I|29|)l and increase the 
threashold value in eq. (|20|l . thus, narrowing the range of 
string parameters relevant for the reionization. If a siz- 
able fraction of matter was ionized by Population III stars 
in the wakes it could lead to a non-monotonic reionization 
similar to the two-step reionization models of 0, |45| . 

In conclusion, our analysis has shown that cosmic 
strings can play a role in early reionization, provided that 
their mass per unit length is within the range given by 
inequalities 120|) and (|29|l . This range is consistent with, 
and may even be somewhat favored by the current data. 
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